Sub-barrier capture with quantum diffusion approach: actinide-based reactions 
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With the quantum diffusion approach the behavior of capture cross sections and mean-square 
angular momenta of captured systems are revealed in the reactions with deformed nuclei at sub- 
barrier energies. The calculated results are in a good agreement with existing experimental data. 
With decreasing bombarding energy under the barrier the external turning point of the nucleus- 
nucleus potential leaves the region of short-range nuclear interaction and action of friction. Because 
of this change of the regime of interaction, an unexpected enhancement of the capture cross section is 
expected at bombarding energies far below the Coulomb barrier. This effect is shown its worth in the 
dependence of mean-square angular momentum of captured system on the bombarding energy. From 
the comparison of calculated and experimental capture cross sections, the importance of quasifission 
near the entrance channel is shown for the actinide-based reactions leading to superheavy nuclei. 
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I. INTRODUCTION 

The measurement of excitation functions down to the 
extreme sub-barrier energy region is important for study- 
ing the nucleus-nucleus interaction as well as the cou- 
pling of relative motion with other degrees of freedom, 
and very little data exist on the fusion, fission and 
ca ptu re cross sections at extreme sub-barrier energies 
[l|-Fl7| . The experimental data obtained arc of inter- 
est for solving astrophysical problems related to nuclear 
synthesis. Indications for an enhancement of the S- 
factor, S = E c , m a exp(2nr]) [HIEII], where rj(E c . m ,) = 
Z\Z2e 2 ' \J /i/(2fi, 2 _B c m .) is the Sommerfeld parameter, at 
energies £c.m. below the Coulomb barrier have been 
found in Refs. 0, E3, EH- However, its origin is still 
under discussion. 

To clarify the behavior of capture cross sections at sub- 
barrier energies, a further development of the theoreti- 
cal methods is required |2(j. The conventional coupled- 
channel approach with realistic set of parameters is not 
able to describe the capture cross sections either below or 
above the Coulomb barrier |13|. The use of a quite shal- 
low nucleus-nucleus potential [2 11 ] with an adjusted repul- 
sive core considerably improves the agreement between 
the calculated and experimental data. Besides the cou- 
pling with collective excitations, the dissipation, which 
is simulated by an imaginary potential in Ref. [2l| or by 
damping in each channel in Ref. (22j , seems to be impor- 
tant. 

The quantum diffusion approach based on the quan- 
tum master-equation for the reduced density matrix has 
been suggested in Ref. |23[. This model takes into con- 
sideration the fluctuation and dissipation effects in colli- 
sions of heavy ions which model the coupling with various 
channels. As demonstrated in Ref. [2J], this approach 
is successful for describing the capture cross sections at 



energies near and below the Coulomb barrier for inter- 
acting spherical nuclei. An unexpected enhancement of 
the capture cross section at bombarding energies far be- 
low the Coulomb barrier has been predicted in [23| . This 
effect is related to the switching off of the nuclear in- 
teraction at the external turning point r ex (Fig. 1). If 
the colliding nuclei approach the distance Rint between 
their centers, the nuclear forces start to act in addition to 
the Coulomb interaction. Thus, at R < Ri n t the relative 
motion is coupled strongly with other degrees of freedom. 
At R > Ri n t the relative motion is almost independent of 
the internal degrees of freedom. Depending on whether 
the value of r ex is larger or smaller than the interaction 
radius Ri n t, the impact of coupling with other degrees of 
freedom upon the barrier passage seems to be different. 

In the present paper we apply the approach of 
Ref. [23} to the description of the capture process 
of deformed nuclei in a wide energy interval includ- 
ing the extreme sub-barrier region. The used for- 
malism is presented in Sect. II. The results of our 
calculations for the reactions 16 0, 19 F, 32 S, 48 Ca+ 232 Th, 
4 Hc, 16 O, 20 Ne, 30 Si, 36 S, 48 Ca+ 238 U, 36 S, 48 Ca, 50 Ti+ 244 Pu, 
48 Ca+ 246 < 248 Cm, and 36 S+ 248 Cm arc discussed in 
Sect. III. The conclusions are given in Sect. IV. 



II. MODEL 

A. The nucleus- nucleus potential 

The potential describing the interaction of two nuclei 
can be written in the form 24 



V(R,Zi,Ai,9i,J) = Vo(R,Zi,Ai,0i, J) 

+ V N (R,Zi,Ai,6i,J) + 
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FIG. 1: (Upper part) The nucleus-nucleus potentials calcu- 
lated at J = (solid curve), 30 (dashed curve), 50 (dotted 
curve), and 65 (dash-dotted curve) for the 16 + 238 U re- 
action. The interacting nuclei are assumed to be spherical in 
the calculation. (Lower part) The position Rb of the Coulomb 
barrier, radius of interaction Ri„t, and external and internal 
turning points for some values of E c . m . are indicated at the 
nucleus-nucleus potential for the same reaction at J—Q. 



the nuclear radius parameter ro = 1.15 fm (for the nuclei with 
At > 16) and the diffuseness parameter a depending on the 
charge and mass numbers of the nucleus [2J]. We use a= 
0.53 fm for the lighter nuclei 16 and 19 F, a= 0.55 fm for 
the intermediate nuclei ( 20 Ne, 26 Mg, 30 Si, 32 > 34 . 36 S, 40 - 48 Ca, 
50 Ti), and a= 0.56 fm for the actinides. For the 4 He nucleus 
r =1.02 fm and a=0.48 fm. 

The Coulomb interaction of two deformed nuclei has the 
following form: 
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where P 2 (cos#i) is the Legendre polynomial. 

In Fig. 1 there is shown the nucleus-nucleus potential V for 
the 16 + 238 TJ reaction (for simplicity, 238 U is assumed to 
be spherical) which has a pocket. With increasing centrifugal 
part of the potential the pocket depth becomes smaller, while 
the position of the pocket minimum moves towards the barrier 

at the position of the Coulomb barrier R — Rb ~ Ri + R2 + 2 

1/3 

fm, where Ri = 1.15 A/ are the radii of colliding nuclei. 
This pocket is washed out at large angular momenta J > 65. 
Thus, only a limited part of angular momenta contributes to 
the capture process. 

For the reactions 36 S + 238 U and 16 O + 238 U (Fig. 2), the 
dependence of the potential energy on the orientation of the 
prolate deformed nucleus 238 U is shown. The lowest Coulomb 
barriers are associated with collisions of the projectile nucleus 
with the tips of the target nucleus, while the highest barriers 
correspond to collisions with the sides of the target nucleus. 
The difference of the Coulomb barriers for the sphere-pole and 
sphere-side orientations is about 16 MeV (8 MeV) for the 36 S 
+ 238 U ( 16 + 238 U) system. 



where Vn, Vc, and the last summand stand for the nuclear, 
the Coulomb, and the centrifugal potentials, respectively. The 
nuclei are proposed to be spherical or deformed. The potential 
depends on the distance R between the center of mass of two 
interacting nuclei, mass Ai and charge Zi of nuclei (i — 1,2), 
the orientation angles 6i of the deformed (with the quadrupolc 
deformation parameters f3i) nuclei and the angular momen- 
tum J. The static quadrupole deformation parameters are 
taken from Ref. [25| for the even-even deformed nuclei. For 
the nuclear part of the nucleus-nucleus potential, we use the 
double-folding formalism, in the form 

V N = J pi(ri)p 2 (R-r 2 )JXri-r 2 )dridr 2 , (2) 

where F(n - r 2 ) = C [F in ^Ill + Fcx( i _ ^)]5(n - r 2 ) 
is the density-dependent effective nucleon-nucleon interac- 
tion and p (r) = pi(r) + p 2 (R - r), F in , ex = / in>ex + 

fin, e J(%l+z\](N 2 2 +z 2 2 ] - Here > Pi( r and N * are the nucleon 
densities and neutron numbers of the light and the heavy nu- 
clei of the dinuclear system, respectively. Our calculations 
are performed with the following set of parameters: Co = 300 
MeV fm 3 , / in = 0.09, / cx = -2.59, f' in = 0.42, = 0.54 and 
poo = 0.17 fm~ 3 [24| . The densities of the nuclei are taken 
in the two-parameter symmetrized Woods-Saxon form with 



B. Capture cross section 

The capture cross section is a sum of partial capture cross 
sections 

^cop(^c.in.) — ^ ^cap(-^c.m., J*) — 
J 

/.tt/2 

= ttA 2 V(2J+1) / d0isin(0i) 
j Jo 

x / d0 2 sin(0 2 )Pcap(£; c . m . I J,0i I 6>2), (4) 
Jo 

where X 2 = ft 2 / (2p_E c . m .) is the reduced de Broglie wave- 
length, p = moAxA^i '(Ai + A2) is the reduced mass (mo is 
the nucleon mass), and the summation is over the possible 
values of angular momentum J at a given bombarding en- 
ergy i?c.m.. Knowing the potential of the interacting nuclei 
for each orientation, one can obtain the partial capture prob- 
ability P cap which is defined by the passing probability of the 
potential barrier in the relative distance R coordinate at a 
given J. 

The value of P cap is obtained by integrating the propagator 
G from the initial state (Ro, Po) at time t = to the final state 
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FIG. 2: (Upper part) The nucleus-nucleus potentials calculated at J—0 for the reactions 36 s_|_ 238 u and 16 + 238 U. (Lower part) 
The dependence of the capture cross section for nuclei colliding with a fixed orientation on E c .m. — y 6 ° rlen * where y° rlent j s the 
height of the Coulomb barrier for certain orientations. The results of calculations for the sphere-sphere (the interacting nuclei 
are spherical), sphere-pole and sphere-side configurations are shown by solid, dashed and dotted lines, respectively. The static 
quadrupole deformation parameters are: /3 2 ( 238 U)=0.286 and /3i( 16 O)=/3i( 36 S)=0. 



(R, P) at time t (P is a momentum): 

/>00 

dR / dP G{R,P,t\Ro,Po,0) 



P cap = lim 

t—¥OQ 



= lim — erfc 



-7-jn + R(t) 



The second line in ([5]) is obtained by using the propa- 
gator G = 7T- 1 |detS- 1 | 1 - /2 e X p(-q T S~ 1 ) (q T = [q R ,qp}, 
q R (t) = R — R(t), q P (t) = P - P(t ), R(t=Q) = Ro, 
P(t = 0) = P , E hk ,(t) = 2q k (t}q k ,(t), S kV (t = 0) = 0, 
k, k' — R, P) calculated in Ref. [2g| for an inverted oscillator 
which approximates the nucleus-nucleus potential V in the 
variable R. The frequency u> of this oscillator with an inter- 
nal turning point r ln is defined from the condition of equality 
of the classical actions of approximated and realistic potential 
barriers of the same hight at given J. It should be noted that 
the passage through the Coulomb barrier approximated by a 
parabola has been previously studied in Refs. [UHH- This 
approximation is well justified for the reactions and energy 
range, which are here considered. Finally, one can find the 



expression for the capture probability: 
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(5) where 7 is the internal-excitation width, ujq 



>{1 



HX'f/[iJ l (si +7)(s2 +7)]} is the renormalized frequency in the 
Markovian limit, the value of A is related to the strength of 
linear coupling in coordinates between collective and internal 
subsystems. The are the real roots (sj > > S2 > S3) of 
the following equation: 



(s + 7)(s 2 — L) ) + hXjs/fj, = 0. 



(7) 



The details of the used formalism are presented in [23j]. We 
have to mention that most of the quantum-mechanical, dis- 
sipative effects and non-Markovian effects accompanying the 
passage through the potential barrier are taken into consid- 
eration in our formalism [23l . [3l| . For example, the non- 
Markovian effects appear in the calculations through the 
internal-excitation width 7. 

As shown in [23|], the nuclear forces start to play a role at 
Rint = Rb + 1-1 fm where the nucleon density of colliding 
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FIG. 3: The calculated value {Vb) averaged over the orien- 
tations of the heavy deformed nucleus versus P c .m. for the 
36 g + 238tj reaction The values of barriers V b orient (sphere- 
pole) for the sphere-pole configuration, Vb = Vb (sphere- 
sphere) = y° rlent (sphere-sphere) for the sphere-sphere config- 
uration and V b orlent (sphere-side) for the sphere-side configu- 
ration are indicated by arrows. The static quadrupole defor- 
mation parameters are: /3 2 ( 238 U)=0.286 and /3i( 36 S)=0. 



nuclei approximately reaches 10% of the saturation density. 
If the value of r ox corresponding to the external turning point 
is larger than the interaction radius Rint, we take Po = r cx 
and Po = in Eq. |5]). For r cx < Ri„t, it is naturally to 
start our treatment with Ro = Rint and Po defined by the 
kinetic energy at R — Ro. In this case the friction hinders 
the classical motion to proceed towards smaller values of R. 
If Po = at Po > Rint, the friction almost does not play a 
role in the transition through the barrier. Thus, two regimes 
of interaction at sub-barrier energies differ by the action of 
the nuclear forces and the role of friction at P = r cx . 
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FIG. 4: The calculated capture cross section (solid lines) ver- 
sus E c . m . for the reactions ie O + 232 Th and 4 He + 238 U are 
compared with the available experimental data. The experi- 
mental data in the upper part are taken from Refs. [37| (o pen 
triangles), [l[ (closed triangles), [3^] (open squares), |3^ | 
(closed squares), [H (open stars) and [5j (closed stars). The 
fission cross sections from Refs. [(| and [36|] are shown in the 
lower part by open circles and solid squares, respectively. The 
value of the Coulomb barrier Vb for the spherical nuclei is indi- 
cated by arrow. The dashed curve represents the calculation 
by the Wong's formula ([8|. The static quadrupole deforma- 
tion parameters are: /3 2 ( 238 U)=0.286, /3 2 ( 232 Th)=0.261 and 
0i( 16 O)=/3i( 4 He)=0. 



III. CALCULATED RESULTS 



Besides the parameters related to the nucleus-nucleus po- 
tential, two parameters hr/— 32 MeV and the friction coeffi- 
cient hX — —h(s\ + S2)—2 MeV are used for calculating the 
capture probability in reactions with deformed actinides. The 
value of A is set to obtain this value of hX. The most real- 
istic friction coefficients in the range of hX « 1 — 2 MeV are 
suggested from the study of deep inelastic and fusion reac- 
tions [32]. These values are close to those calculated within 
the mean field approach [33l ]. All calculated results presented 
are obtained with the same set of parameters and are rather 
insensitive to a reasonable variation of them [23l . |2Sl l3l| . 



A. Effect of orientation 

The influence of orientation of the deformed heavy nucleus 
on the capture process in the reactions 36 S + 238 U and 16 O 
_l_ 238tj j s s t u died in Fig. 2. We demonstrate that the cap- 
ture cross section a cap at fixed orientation as a function of 
Pc.m. - V b orient , where V h orimt is the Coulomb barrier for this 
orientation, is almost independent of the orientation angle (9 2 ■ 

In Fig. 3 the value of the Coulomb barrier 

<V b > = -^](2J + 1) f ' d9 2 sm(e 2 ) 

& cap (-^c.m.J j J Q 

X Pca.p(E c .m., J,Ol,02)V{Rb, Zi, Ai,0i, J) 
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FIG. 5: The same as in Fig. 4, but for the reactions O 
+ 238 U and 36 S + 238 U. The experimental cross sections are 
taken from Refs. (01 (open triangles), I35l| (closed triangles) , 
[3 (open squares), [3| ( closed squares), [3(| (open stars), 
(closed stars), and [14| (rhombuses). The dashed curve rep- 
resents the calculation by the Wong's formula ((5J. The static 
quadrupole deformation parameters are: /32( 238 U)=0.286 and 
/3i( 16 O)=/M 36 S)=0. 



averaged over all possible orientations of the heavy nucleus 
versus E c , m . is shown for the 36 S + 238 U reaction. With in- 
creasing (decreasing) E c . m . the value of < Vb > approaches 
the value of the Coulomb barrier for the sphere-sphere config- 
uration (for the sphere-pole configuration). The influence of 
deformation on the capture cross section is very weak already 
at bombarding energies about 15 MeV above the Coulomb 
barrier corresponding to spherical nuclei. 



B. Comparison with experimental data and 
predictions 

In Figs. 4-6 the calculated capture cross sections for the re- 
actions ie O, 19 F, 32 S+ 232 Th and 4 He, 16 O, 30 Si, 32 ' 36 S+ 238 U are 
in a rather good agreement with the available experimental 
data [J, d |1S, & E3, USEE HUU. Because of the uncer- 
tainties in the definition of the deformation of the light nu- 
cleus and in the experimental data [HI . l4(il ) in Fig. 6, we 
show the calculated results for the 30 Si+ 238 U reaction with 



FIG. 6: The same as in Fig. 4, but for the reactions ' S 
+ 232 Th (solid line), 32 S + 238 U (dotted line) and 30 Si + 
238 U. The experimental data are taken from Refs . 1151 ] ( 32 S 
_l_ 232 rpj^ solid squares), [ll|] (solid circles) and [4Q] (open 
squares). The static quadrupole deformation parameters 
are: /3 2 ( 238 U)=0.286, /3 2 ( 232 Th)=0.261, /3i ( 32 S)=0.312 and 
/3i( 3O Si)=0.315. For the 30 Si + 238 U reaction, the results of 
calculations with /3i( 30 Si)=0 (the predictions of the mean- 
field and macroscopic-microscopic models) are presented by 
dashed line in the lower part of the figure. 



/3i( 30 Si) from Ref. [H] as well as with /3i( 3o Si)=0 (lower part 
of Fig. 6). Note that /3i( 30 Si)=0 for the ground state were 
predicted within the mean-field and macroscopic-microscopic 
models. 

In Fig. 7 (upper part) we are not able to describe well 
the data of Ref. [| for the 19 F+ 232 Th reaction at E c . m . < 
74 MeV, even by varying the static quadrupole deformation 
parameters /3i of 19 F. However, the deviations of the solid 
curve in the upper part of Fig. 7 from the experimental data 
are within the uncertainty of these data. Note that the value 
of /?i mainly influences the slope of curve at E c .m. < Vb and 
one can extract the ground state deformation of nucleus from 
the experimental capture cross section data. For the 20 Ne + 
238 U reaction, the calculated capture cross sections in Fig. 7 
are consistent with the experimental data [3(| if the latter 
ones are shifted by 2 MeV to lower energies. For the 20 Ne 
nucleus, the experimental quadrupole deformation parameter 
/3i =0.727 related in Ref. [25|] to the first 2 + state seems to 
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FIG. 7: The same as in Fig. 4, but for the reactions F + 
232 Th and 20 Ne + 238 U. The experimental data are taken from 
Refs. [f| (open squares), Q| (closed squares), and [3(| (closed 
circles). The open circles in the lower part are the experi- 
mental data from Ref. shifted by 2 MeV to the left. The 
results of calculations with the static quadrupole deformation 
parameters of 19 F /3i ( 19 F)=0.275 (as in Ref. 0), 0.41, and 
0.55 are shown by the dashed, solid, and dotted lines, respec- 
tively. The other static quadrupole deformation parameters 
are: /3 2 ( 238 U)=0.286, /3 2 ( 232 Th)=0.261 and /3i( 20 Ne)=0.335. 



be unrealistically large and we take f3\ =0.335 as predicted in 
Ref. [4l[. The capture cross sections for the reactions 32 S 
+ 238 U and 36 S + 244 Pu, 248 Cm are shown in Figs. 6 and 8, 
respectively. 

One can see in Figs. 4-8 that there is a sharp fall-off of the 
cross sections just under the Coulomb barrier corresponding 
to undeformed nuclei. With decreasing E c . m . up to about 
8-10 MeV (when the projectile is spherical) and 15-20 MeV 
(when both projectile and target are deformed nuclei) below 
the Coulomb barrier the regime of interaction is changed be- 
cause at the external turning point the colliding nuclei do 
not reach the region of nuclear interaction where the friction 
plays a role. As result, at smaller E c .m. the cross sections fall 
with a smaller rate. With larger values of Rint the change 
of fall rate occurs at smaller E c . m .. However, the uncertainty 
in the definition of Rint is rather small. Therefore, an effect 
of the change of fall rate of sub-barrier capture cross section 
should be in the data if we assume that the friction starts 



FIG. 8: The predicted capture cross sections for the reac- 
tions 36 S + 244 Pu, 248 Cm. The static quadrupole deforma- 
tion parameters are: /3 2 ( 244 Pu)=0.293, /3 2 ( 248 Cm)=0.297 and 
/3i( 36 S)=0. 



to act only when the colliding nuclei approach the barrier. 
Note that at energies of 10-20 MeV below the barrier the ex- 
perimental data have still large uncertainties to make a firm 
experimental conclusion about this effect. The effect seems 
to be more pronounced in collisions of spherical nuclei, where 
the regime of interaction is changed at -E c .m. up to about 3-5 
MeV below the Coulomb barrier [2^|. The collisions of de- 
formed nuclei occur at various mutual orientations affecting 
the value of Rint- 

The well-known Wong formula for the capture cross section 

is 

D2t rn/2 nn/2 

a(E c . m .) = -J* / d6>isfn0i / d6 2 sm0 2 

2-C/c.m. J() JO 

x ln(l + exp[27r(£ c .m. - -E 6 (0i, 2 ))/M)> (8) 

where _Ei,(#i,# 2 ) is value of the Coulomb barrier which de- 
pends on the orientations of the deformed nuclei [44| . As seen 
from Figs. 4 and 5 (dashed lines) the Wong formula ((8J does 
not reproduce the capture cross section at E c . m , < Vb even 
taking into consideration the static quadrupole deformation 
of target-nucleus. 
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FIG. 9: The calculated mean-square angular momenta ver- 
sus E a .m. for the reactions ls O + 232 Th, 238 U are compared 
with experimental data The dashed curve represents 

the calculation by the Eq. (|10[) . The static quadrupole de- 
formation parameters are: /3 2 ( 238 U)=0.286, /3 2 ( 232 Th)=0.261 
and /3i( 16 O)=0. The values of the Coulomb barriers VJ, cor- 
responding to spherical interacting nuclei are indicated by 
arrows. 



FIG. 10: The same as in Fig. 9, but for the indicated reac- 
tions 19 F, 48 Ca + 232 Th. The experimental data are taken 
from Ref. 0]. The results of calculations with quadrupole 
deformation parameters f3\ ( 19 F)=0.275, 0.41 and 0.55 are 
shown by the dashed, solid, and dotted lines, respectively. 
The other static quadrupole deformation parameters are: 
/3 2 ( 232 Th)=0.261 and /3i( 48 Ca)=0. 



The calculated mean-square angular momenta 
(J 2 ) = 



nX 2 J2 j J(J + 1)(2 J +1) 



u cap (Bern.) 
/2 r-TT/2 

d0isin(0i) / dO2sm(0 2 )P C!lp (E c . m .,J,9i,02) 

(9) 

of captured systems versus -E c .m. are presented in Figs. 9-10 
for the reactions mentioned above. At energies below the bar- 
rier the value of (J 2 ) has a minimum. This minimum depends 
on the deformations of nuclei and on the factor Z\ XZ2. For 
the reactions 16 + 232 Th, le O + 238 U, 19 F + 232 Th and 48 Ca 
_l_ 232^^ t nese minima are about 7, 8, 12 and 15 MeV below 
the corresponding Coulomb barriers, respectively. The ex- 
perimental data [43I indicate the presence of the minimum as 
well. On the left-hand side of this minimum the dependence 



of (J ) 



has been found in Refs. [43j] in the extreme sub-barrier region. 
Note that the found behavior of (J 2 ), which is related to the 
change of the regime of interaction between the colliding nu- 
clei, would affect the angular anisotropy of the products of 
fission-like fragments following capture. Indeed, the values 
of (J 2 ) are extracted from data on angular distribution of 
fission- like fragments [l7j |. 

In the Wong model [44| the value of the mean-square an- 
gular momentum is determined as 



(J ) 



irh 2 



t/2 



dOi sin 0i 



ir/2 



d02 sin 02 



-L^ 2 (-exp[27r(£ c . m . - Eb(0i , 6> 2 ))/H) , 1fA 
X - — — — — — — — • (.tUJ 



ln(l + exp[27r(_B c 



Si,(6>i,6»2))/M) 



is rather weak. A similar weak dependence rier 



Here, the Li2^z) is the polylogarithm function. At 
exp[27r(_E c m . — Et)/hu}]) <Sl (much below the Coulomb bar- 

Li 2 (- C xp[2^(E c . m .-E,,)/M) ^ 1 and Qne can obtain 



ln(l + cxp[2ir(B c 



-B 6 )/M) 
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FIG. 11: The calculated values ol the astrophysical S- 
factor with 770 = rj(E c . m . = 14) for the indicated reactions 

16 0+ 232 Th 4 He + 238U q£ ^ Coulomb bar _ 

riers Vj, corresponding to the spherical nuclei are 78.6 and 21.2 
MeV. 



saturation value of the mean-square angular momentum [2 



<J 2 ) = 



nh 2 



(11) 



The agreement between (J 2 ) calculated with Eq. (|10p and 
experimental (J 2 ) is not good. At energies below the barrier 
{J 2 } has no a minimum (see Fig. 9). However, for the consid- 
ered reactions the saturation values of (J 2 ) are close to those 
obtained with our formalism. 



Astrophysical factor, L-factor and barrier 
distribution 



In Figs. 11 and 12 the calculated astrophysical S-factors 

versus E c .m. are shown for the reactions 4 He, 16 0+ 238 U and 
i6 0+ 232 Th The S _ factor has 

a maximum for which there 
are experimental indications in Refs. 0, fiol . l2ll ]. After this 
maximum S-factor slightly decreases with decreasing E c . m . 
and then starts to increase. This effect seems to be more 



FIG. 12: The calculated values of the astrophysical S-factor 
with 770 = T?(-E c .m. = Vb) (middle part), the logarithmic 
derivative L (upper part) and the fusion barrier distribu- 
tion d 2 (£ c .m.cr CQ p)/d_B 2 m . (lower part) for the 16 + 238 U re- 
action. The value of L calculated with the assumption of 
/3i( 16 O)=/3 2 ( 238 U)=0 is shown by a dashed line. The solid 
and dotted lines show the values of d 2 (-E c .m.o" cap ) /dE 2 m cal- 
culated with the increments 0.2 and 1.2 MeV, respectively. 
The closed squares are the experimental data of Ref. [50 ]. 
The value of the Coulomb barrier 14 corresponding to the 
spherical nuclei is 80 MeV. 



pronounced in collisions of spherical nuclei [23|. The same 
behavior has been revealed in Refs. [34| by extracting the 
5-factor from the experimental data. 

In Fig. 12, the so-called logarithmic derivative, 
i(-E c .m.) = d(ln(E c . m a cap ))/dE c . m ., and the barrier 
distribution d 2 {E c . ul .a cav )/dE 2 . ra . are presented for the 
16 + 238 U reaction. The logarithmic derivative strongly 
increases below the barrier and then has a maximum at 
E c . m . ~ V 6 orle " t (sphere-pole)-3 MeV (at E c . m . w 14-3 
MeV for the case of spherical nuclei). The maximum of L 
corresponds to the minimum of the S-factor. 

The barrier distributions calculated with an energy in- 
crement 0.2 MeV have only one maximum at -E c . m . « 
V b orient (sphere-sphere) = V b as in the experiment With 
increasing increment the barrier distribution is shifted to 
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FIG. 13: The same as in Fig. 4, but for the 48 Ca + 2 3 2 Th, 238 U 
reactions. The excitation energies Eq N of the correspond- 
ing nuclei are indicated. The experimental data are taken 
from Refs. (marked by squares) and [47| (marked by 
circles). The static quadrupole deformation parameters are: 
/3 2 ( 238 U)=0.286, /3 2 ( 232 Th)=0.261 and /3i( 48 Ca)=0. 



FIG. 14: The same as in Fig. 13, but for the reactions Ca 
+ 246,248 Cm. The experimental data are from Refs. 01 
(squares) and [4S|] (circles). The static quadrupole defor- 
mation parameters are: /3 2 ( 246 Cm)=0.298, /3 2 ( 248 Cm)=0.297 
and /3i( 48 Ca)=0. 



lower energies. Assuming a spherical target nucleus in the 
calculations, we obtain a more narrow barrier distribution 
(see Fig. 12). 



D. Capture cross sections in reactions with large 
fraction of quasifission 

In the case of large values of Z\ x Z 2 the quasifission pro- 
cess competes with complete fusion at energies near barrier 
and can lead to a large hindrance for fusion, thus ruling the 
probability for producing superheavy elements in the actinide- 
based reactions [45l . [4r3 ] . Since the sum of the fusion cross 
section <r/ us , and the quasifission cross section a q f gives the 
capture cross section, 

f cap — G fus &qf , 



and Uf us <§; cT q f in the actinide-based reactions 48 Ca + 
232 Th, 238 U, 244 Pu, 246 ' 248 Cm and 50 Ti + 244 Pu [H, we have 



" cap O qf . 

In a wide mass-range near the entrance channel, the quasi- 
fission events overlap with the products of deep-inelastic col- 
lisions and can not be firmly distinguished. Because of this 
the mass region near the entrance channel is taken out in 
the experimental analyses of Refs. ptsl . |49| . Thus, by com- 
paring the calculated and experimental capture cross sections 
one can study the importance of quasifission near the entrance 
channel for the actinide-based reactions leading to superheavy 
nuclei. 

The capture cross sections for the quasifission reactions [47l — 
[4S| are shown in Figs. 13-15. One can observe a large devia- 
tions of the experimental data of Refs. from the the 
calculated results. The possible reason is an underestimation 
of the quasifission yields measured in these reactions. Thus, 
the quasifission yields near the entrance channel are impor- 
tant. Note that there are the experimental uncertainties in 
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FIG. 15: The same as in Fig. 13, but for the indicated 
48 Ca, 50 Ti + 244 Pu reactions. The experimental data are 
from Refs. [49l ] (squares) and [48|] (circles). The static 
quadrupole deformation parameters are: /?2( 244 Pu)=0.293, 
and / 9i( 48 Ca)= J 8i( B0 Ti)=0. 




20 30 40 50 

E* CN (MeV) 



the bombarding energies. 

One can see in Fig. 16 that the experimental and the theo- 
retical cross sections become closer with increasing bombard- 
ing energy. This means that with increasing bombarding en- 
ergy the quasifission yields near the entrance channel mass- 
region decrease with respect to the quasifission yields in other 
mass-regions. The quasifission yields near the entrance chan- 
nel increase with Z\ x Zi. 



IV. SUMMARY 

The quantum diffusion approach is applied to study the 
capture process in the reactions with deformed nuclei at sub- 
barrier energies. The available experimental data at ener- 
gies above and below the Coulomb barrier are well described, 
showing that the static quadrupole deformations of the inter- 
acting nuclei are the main reasons for the capture cross section 
enhancement at sub-barrier energies. Since the deformations 
of the interacting nuclei mainly influence the slope of curve 
at E c . m , < Vt and one can extract the ground state defor- 
mation of projectile or target from the experimental capture 
cross section data. 

Due to a change of the regime of interaction (the turning- 
off of the nuclear forces and friction) at sub-barrier energies, 
the curve related to the capture cross section as a function of 
bombarding energy has smaller slope E c . m . — V& < - 5 MeV. 
This change is also reflected in the functions (J 2 ), L(E c . ln .), 
and S(E c .m.). The mean-square angular momentum of cap- 
tured system versus E c . m . has a minimum and then satu- 
rates at sub-barrier energies. This behavior of (J 2 ) would 
increase the expected anisotropy of the angular distribution 
of the products of fission and quasifission following capture. 
The astrophysical factor has a maximum and a minimum at 
energies below the barrier. The maximum of L-factor corre- 
sponds to the minimum of the S-factor. One can suggest the 
experiments to check these predictions. 

The importance of quasifission near the entrance channel is 
shown for the actinide-based reactions leading to superheavy 
nuclei. 



FIG. 16: The ratio of theoretical and experimental capture 
cross sections versus the excitation energy E c . m . of the com- 
pound nucleus for the reactions 4S Ca+ 238 U (closed stars), 
48 Ca+ 244 Pu (closed triangles), 48 Ca+ 246 Cm (closed squares), 
48 Ca+ 248 Cm (closed circles), and 50 Ti+ 244 Pu (closed rhom- 
buses). 
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